Precision Feshbach spectroscopy of ultracold Cs2 
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We have observed and located more than 60 magnetic field-induced Feshbach resonances in ultra- 
cold collisions of ground-state ^'^''Cs atoms. These resonances are associated with molecular states 
with up to four units of rotational angular momentum, and are detected through variations in the 
elastic, inelastic, and radiative collision cross sections. These observations allow us to greatly im- 
prove upon the interaction potentials between two cesium atoms and to reproduce the positions of 
most resonances to accuracies better than 0.5%. Based on the relevant coupling scheme between 
the electron spin, nuclear spin, and orbital angular momenta of the nuclei, quantum numbers and 
energy structure of the molecular states beneath the dissociation continuum are revealed. Finally, 
we predict the relevant collision properties for cesium Bose-Einstein condensation experiments. 
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I. INTRODUCTION 

The collision properties of an ultracold and dilute atom 
gas are strongly influenced by the long-range interactions 
between two atoms. When the interaction potential sup- 
ports a weakly bound state near the scattering energy, 
the atomic collision properties can be resonantly altered, 
a situation referred to as a Feshbach resonance P, 
In many cold-atom systems, magnetically tunable Fesh- 
bach resonances have been discovered, and have led to 
ground-breaking observations including the implosion of 
a Bose-Einstein condensate (BEC) the coherent cou- 
pling between an atomic BEC and molecules Q , the cre- 
ation of bright solitonsja, and recently the creation of 
ultracold molecules |^ UlEl a-nd of a molecular BEC 9] . 

The collisional properties of ultracold cesium atoms 
have intrigued experimentalists and theorists because of 
their large clock shifts |1^J , enormous collision cross sec- 
tions 0, and extreme difficulty to reach BEC [l^ . 
These anomalies in the atom-atom scattering can be ex- 
plained by the coupling of the scattering continuum to 
molecular states. While these states cannot be accessed 
by conventional spectroscopy they may be tuned 
into resonance with the scattering continuum and induce 
Feshbach resonances. Detection of multiple Feshbach res- 
onances, or Feshbach spectroscopy, permits a precise de- 
termination of the long-range interaction parameters, as 
well as the molecular structure near threshold. With this 
information, the cold collision anomalies can be resolved 



and the clock shifts 14], collision cross sections and scat- 
tering lengths can be accurately calculated 0. 

In this work, we report the observation of more than 
60 Feshbach resonances of cesium atoms in 10 different 
incident channels. In particular, we employ a radiative 
detection scheme to resolve narrow resonances 0, IT^ 
whose locations allow us to significantly improve our de- 
termination of Cs interaction parameters over our previ- 
ous work With these parameters, the molecular 
energy structure near threshold, as well as s-wave scat- 
tering lengths and collision properties can be precisely 
determined. 

We organize the paper as follows. First, we outline 
the experimental setup and procedures, and the general 
measurement methods in SecE] We present and discuss 
the results from inelastic Feshbach spectroscopy, elastic 
Feshbach spectroscopy and radiative spectroscopy in Sec- 
tions mil IIVI and respectively. In Section IVII we dis- 
cuss the Hamiltonian for two ground-state cesium atoms 
and approximate quantum numbers in the system. The 
numerical procedures for calculating scattering proper- 
ties and bound state energies are presented in Sec. IVIII 
Section IVIIII analyses all observed Feshbach resonances 
and assigns quantum numbers to each resonance. 



II. EXPERIMENT 



A. Feshbach spectroscopy 
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Feshbach resonance in binary atomic collisions is illus- 
trated in Fig. ^ Interacting atom pairs in the scatter- 
ing continuum, or scattering channel, couple to a discrete 
bound state supported by a closed channel with higher in- 
ternal energy. This coupling resonantly alters the outgo- 
ing scattering amplitude in the scattering channel (elastic 
collision resonance) and in the lower-lying open channels 
(inelastic collision resonance). 
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interatomic separation 

FIG. 1: Illustration of radiative, elastic, and inelastic Fes- 
hbach resonance. The molecular state is supported by the 
closed channel. 



In ultracold atom experiments, Feshbach resonances 
can be observed in elastic and inelastic collision rates 
0, E| and in the molecular population, probed by pho- 
toassociation or photodissociation 0| as shown in 
Fig. [2 

In most cases the magnetic moments of the molecu- 
lar state and the scattering continuum are unequal, such 
that the energy difference between the molecular state 
and the colliding atom pair can be tuned by means of 
an external magnetic field. The molecular bound-state 
energy structure near dissociation threshold is then re- 
flected in the magnetic-field dependence of the atom- 
atom scattering. A measurement of the ultracold colli- 
sion properties as a function of magnetic field (Feshbach 
spectrum) in combination with theoretical modeling thus 
unveils the underlying molecular spectrum. In this work, 
we are able to associate each resonance with its quantum 
numbers of both the incoming scattering channel and the 
molecular state, where it is noteworthy that the s-, p-, d- 
. . . partial wave of the incoming channel and of the res- 
onance need not be the same. A collection of Feshbach 
spectra for colliding atoms in several different quantum 
states then provides the information necessary to accu- 
rately determine the long-range interaction parameters 
and the molecular energy structure near the dissociation 
continuum. 



B. Preparation of ultracold high-density samples 

In order to measure Feshbach spectra with good signal- 
to-noise ratio and high resolution, and to simplify the 
theoretical analysis, it is favorable to use high-density 
samples of ultracold atoms confined in a trap that pro- 
duces negligible tensor light shifts of the atomic and 
molecular energy levels [23 ■ We use a far-detuned, lin- 
early polarized, one-dimensional (ID) optical lattice trap 
to confine the atoms, in conjunction with an optical cool- 
ing method that c an p roduce high-density samples at fiK 
temperatures 
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FIG. 2: Apparatus for Feshbach spectroscopy of ^^''Cs. 



Fig.Elshows schematically the experimental setup. Ul- 
tracold samples of cesium atoms are prepared in the up- 
per part of a compact ultrahigh vacuum chamber. A 
pair of magnetic field coils, when operated in the anti- 
Helmholtz configuration with opposite currents, provides 
the spherical quadrupole field for a magneto-optical trap 
(MOT), and, in the Helmholtz configuration, produces 
a homogeneous magnetic field up to 25 mT. A linearly 
polarized, vertically propagating, retroreflected Nd:YAG 
laser beam provides the far-detuned one-dimensional lat- 
tice dipole trap at 1064 nm j2l|. The number of trapped 
atoms is inferred from the fluorescence emitted by the 
cloud when it is illuminated with resonant light on the 
Fg = 4: —> Ff, = 5 hyperfine component of the D2 line 
near 852 nm. Cloud temperatures in the vertical and hor- 
izontal directions are determined by releasing the atoms 
and performing time-of-flight imaging 12 cm below onto 
a photodiode and a linear CCD array, respectively. Trap 
vibration frequencies are measured by parametric excita- 
tion |2,l| . The atomic density is derived from the mea- 
sured atom number, temperature and trap vibration fre- 
quencies, as detailed in Ref. 01 ■ 

The ultracold sample is prepared by first collecting 
5 X 10^ atoms in a vapor-cell magneto-optical trap (MOT) 
in 500 ms. Superimposed with the MOT is the YAG 
ID optical lattice dipole trap that at a power of 8 W 
and a beam waist of 260 /^m provides a trap depth of 
U/h — 1.6 MHz, and axial and radial vibration frequen- 
cies of oJa/'^Tr = 50 kHz and utr/^n — 80 Hz, respectively. 
The preparation of the high-density sample in the YAG 
dipole trap is accomplished by means of two phases of 
Raman-sideband coohng (RSC) with a 3D near-detuned 
optical lattice [2^. We extinguish the MOT light and 
apply the first RSC for 10 ms in a small bias field of 5 
/iT, which cools the atoms to temperatures below 1 /iK, 
and predominantly polarizes them into the lowest-energy 
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magnetic sublevel |6S'i/2,F = 3,mp = 3). Here F is the 
total angular momentum of cesium atoms in the 6S'i/2 
ground state and mp is the angular momentum projec- 
tion along the magnetic field direction. When the 3D 
near-detuned optical lattice is extinguished in 1 ms, the 
release of the atoms into the II? YAG lattice trap is adi- 
abatic only in the vertical direction, and we observe a ra- 
dial oscillation as the atoms slide down the trap potential. 
To remove some of the excess potential energy, we wait 
for 4 ms until the atoms have the greatest kinetic energy 
and highest density, and perform a second phase of RSC 
p^. Due to the high density near 10"'^^ cm~'^, the second 
cooling requires a weaker optical pumping intensity and 
longer cooling time of 15 ms. After rethermalization for 
200 ms in the YAG lattice trap, 1 x 10* atoms are pre- 
pared at a temperature of 3 fiK to 5 fiK with a vertical 
(horizontal) rms radius of cr^ = 580 /im {ax = 30 /im). 
This means that near the central region, each site in the 
ID YAG lattice trap contains 4 x 10* atoms at a mean 
density of approximately 1 x lO^^cm"'^. The two phases 
of RSC also allow us to adjust almost independently the 
atomic density and temperature by changing the detun- 
ing of the optical pumping beam during the first and the 
second RSC phase, respectively. 

The preparation of pure samples in desired target 
states is crucial for cold collision experiments because 
different internal states can have drastically different col- 
lision properties. To improve the atomic polarization 
over that achieved by the RSC alone, we apply an ad- 
ditional optical pumping pulse for 3 ms during the re- 
lease from the second RSC phase. This has the advan- 
tage that atom-atom collisions are still suppressed due 
to the SD confinement in the near-detuned optical lat- 
tice, which prevents radiative collision loss. Using mi- 
crowave spectroscopy, we optimize the optical pumping 
process at a bias field of 14 /iT, and prepare up to 98% of 
the atoms in one of the stretched states |-F, tof = 
and the remaining atoms in the neighboring \F,mp = 
±{F — 1)) state, where F—3 or 4. Fcshbach resonances 
in four scattering states or channels can be thus stud- 
ied: (3,3) + (3,3), (3,-3) + (3,-3), (4,4) + (4,4) and 
(4, —4) + (4, —4), where {Fi,mpi) + {F2,mp2) indicates 
the collision of one atom in the |-Fi, mpi) state with one 
in the \F2,mp2) state. 

We also prepare mixed samples containing atoms in 
two different internal states by either detuning the optical 
pumping beam, or by applying an additional microwave 
pulse to transfer part of the population. Typically, we 
prepare 90% of the population in a stretched state |a), 
and 10% in another state |6). After fully characterizing 
the collision properties of the (a) 4- (a) channel, the mixed 
collisions (a) -t- (6) can be monitored by either identifying 
resonances growing with the population in or by se- 
lectively detecting only the population in The large 
population ratio of the two states ensures that (6) -I- (6) 
collision processes remain insignificant. Mixed scattering 
channels investigated in this work include (3, 3) -I- (3, 2), 
(3, 3) + (4, 4), (3, 3) + (4, 3), (3, 3) + (4, 2),(3, -3) + (3, -2) 



and (4,-4)-f (4, -3). 

The ability to study all these different scattering chan- 
nels is crucial to provide information on the molecular 
bound states near the (Fi = 3) + {F2 = 3), {Fi = 
3) + {F2 = 4), and (Fi = 4) + (F2 = 4) hyperfine asymp- 
totes. These span roughly h 18 GHz in binding energy. 
As was shown in Ref. (15|i the molecules have different 
characteristics near the three dissociation limits, which 
particularly helps determine the strength of the van der 
Waals interaction. 

After the atoms have thermalized in the dipole trap, 
a uniform magnetic field up to 25 mT is applied. In or- 
der to preserve the atomic polarization during the field 
ramp, we first increase the field from 14 /xT to 200 /j,T 
in 200 ms, and then to an arbitrary field value in an- 
other 100 ms. The magnetic field experienced by the 
atoms is calibrated to an accuracy of SB < O.l/iT at low 
field and SB/B < 10~^ at high field from the Zeeman 
splitting between magnetic sublevels, as measured with 
microwave spectroscopy . Slowly changing stray fields 
of typically 50 f/T and remnant fields from the magne- 
tized vacuum chamber up to 2 /iT are carefully cancelled 
with six independent bias coils in three orthogonal di- 
rections to an accuracy of 0.1 fiT for the field, and 2 
/iT/cm for the field gradient. An effective magnetic field 
due to residual circular polarization of the YAG trapping 
beam |20j is monitored with microwave spectroscopy, and 
reduced below 0.05 /xT by linearizing the beam polariza- 
tion. To suppress the effects of field inhomogeneity and 
atomic density variation between the various YAG ID 
lattice sites, we perform measurements only on the cen- 
ter portion of the cloud within ±0.4(72, which contains 
900 lattice sites with a mean atomic density variation of 
< 10%. 



III. INELASTIC FESHBACH SPECTROSCOPY 

Inelastic collisions occur when the initial scattering 
state couples to open channels with lower internal en- 
ergy, see Fig. ^ Due to energy conservation, the internal 
energy difference is converted into relative kinetic energy 
of the atom pair, and is either on the order of the ground 
state hyperfine splitting or the Zeeman energy. In both 
cases, the energy release is generally much larger than 
our trap depth and therefore results in a loss of the col- 
liding atom pair from the trap. Near a Feshbach reso- 
nance, inelastic loss is either enhanced or suppressed due 
to constructive or destructive interference between the 
off-resonant scattering amplitude with the on-resonant 
amplitude. In this work, we observe mostly enhanced 
inelastic collision processes, and only one prominent sup- 
pression at 7^6 mT in the (3, 3) -I- (4, 2) incident scatter- 
ing channel (tJ- 

Experimentally, the inelastic rates are determined by 
fitting the mean atom density, n, to dn/dt = —Ln — Kn^, 
where L is an one-body loss rate, and K is the desired 
collision rate coefficient. If we assume that the atoms 
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FIG. 3: Resonant and ofF-resonant time evolution of the 
atomic density. Atoms are prepared in |_F = 4, mp — —4) 
at T = 5.3 fiK and n = 2.1 x 10^^ cm~^ at B = 14 fiT (solid 
circles) and B — 20.53 mT (open circles), where a strong 
Feshbach resonance is located. See also Table |T| 



A. (4,4) + (4,4) inelastic collisions 

The stretched |4,4) state was considered a promising 
candidate for reaching Bose-Einstein condensation in a 
magnetic trap before the cesium collision properties were 
revealed. Large inelastic losses that grow as the temper- 
ature is reduced were discovered that prevent the con- 
densation of cesium in this state jl^, 123 • 

We find no resonance structure for a sample polarized 
in the |4,4) state at 5/iK, but a loss rate coefficient in- 
creasing slowly from 2x 10~^^ cm^s~^ to 3x 10~^^ cm^s~^ 
over the range of i? = to 25 mT. The absence of res- 
onances in this channel is expected since the scattering 
channel (4, 4) -I- (4, 4) has the highest hyperfine and Zee- 
man energy of all ground-state hyperfine levels, and con- 
sequently Feshbach resonances cannot occur. 



B. (4,-4) + (4,-4) and (4,-4) + (4,-3) inelastic 
collisions 



only reside in the harmonic region of the YAG dipole trap 
potential and that the cloud is in thermal equilibrium 
and at constant temperature T during the experiment, 
then the mean atom density is n = N {muj'^ / iTrkBT)^^^ 
where N is the atom number in one lattice site, uj — 

1/32/3. 

LJz i-^r is the mean vibration frequency, m is the cesium 
atomic mass, and ks is the Boltzmann constant. The 
above differential equation can be solved analytically and 
directly relates the measured trap loss to the thermally 
averaged loss coefficient K. The one-body loss rate is 
the same for all measurements. A typical measurement 
is shown in Fig. |31with a holding time up to 5 s. 

The loss rates as a function of external magnetic field 
are measured by observing the atom loss within a holding 
time between 30 ms to 300 ms, such that the maximum 
collision loss is less than 30%, and the atomic temper- 
ature varies by less than 10%. Slow variations in the 
initial atom number are monitored after every 10 mea- 
surements, and corrected for. On the other hand, to 
observe and locate weak resonances, we allow the atoms 
to interact for a longer time up to 500 ms in order to 
obtain better signal to noise. In this case, we ignore the 
temperature evolution and report only the fractional loss 
of atoms. 

Collisions between atom pairs in two different internal 
states, such as (3, —3) -|- (3, —2), differ from those in the 
same state in that both the even and odd partial waves 
can be scattered. In the (3, —3) -I- (3, —2) channel, for in- 
stance, we observe rich Feshbach spectra containing both 
s- and p-wave resonances. 

Positions of the inelastic Feshbach resonances are tab- 
ulated in Tables HI and ITU 



We observe two narrow inelastic Feshbach resonances 
for collisions between two |4, —4) atoms: a weak one at 
10.59 mT and a strong one at 20.50 mT. In addition we 
have discovered an inelastic resonance in the (4, —4) -|- 
(4, —3) channel at 20.66 mT. The resonance field values 
have a ratio of 1 : 2 : 2 to within 3%. This ratio is not 
a coincidence, but a result of the molecular bound state 
structure below the continuum. See Table All three 
resonances are identified as originating from bound states 
with identical binding energy at zero magnetic field. 



C. (3,-3) + (3,-3) and (3,-3) + (3,-2) inelastic 
collisions 

The |3, —3) state can be magnetically trapped and sev- 
eral attempts to reach EEC in this state were thwarted 
by the large collision loss and negative scattering length 
at low magnetic field |23|. In this work, we observe mul- 
tiple inelastic resonances for collisions between two such 
atoms, as shown in Fig.^ One resonance at 2.18 mT is 
not identified, see Table ITU 

In (3, —3) + (3, —2) collisions, multiple s- and p- wave 
resonances are found and identified by varying the pop- 
ulation in the |3, —2) state, shown in Fig. |5| An alterna- 
tive method to identify mixed state resonances is based 
on detecting the population in |3, —2) by microwave tran- 
sitions on samples with 90% of the population in |3, —3) 
and ~ 10% in |3,— 2). In this case, we observe enhanced 
loss of the |3, —2) population near inelastic collision res- 
onances in the (3, —3) -f (3, —2) channel, and suppressed 
loss of the |3, —2) population near inelastic collision res- 
onances in the (3, —3) -I- (3, —3) channel, the latter being 
due to the reduction in density of |3,— 3) atoms. This 
provides a clear distinction between the two collision pro- 
cesses. 
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FIG. 4; Binary loss coefficient in a gas of|F = 3,mF = — 3) 
atoms as a function of external magnetic field. The initial 
mean atomic density is n = 5 x lO^'^ cm~^ at a temperature 
of ~ 5 /iK. Populations in the \F — 3,mF ~ —3), \F — 
3, rriF = —2) and all other states are 95%, ~ 5% and < 1%, 
respectively. Magnetic field resolution is lO/itT. 




magnetic field B [ml] 

FIG. 5: Atom loss as a function of magnetic field for an initial 
mean density n = 5 x lO^'^cm"'^, temperature of ~ 5 fiK and 
interaction time of 300 ms. Populations in the \F = 3,mF ~ 
—3), \F = 3,mF = —2) and all other states are 95%, ~ 5% 
and < 1% for the lower curve, 85%, ~ 10% and ~ 5% for the 
middle curve and ~ 70%, ~ 20% and ~ 10% for the upper 
curve, respectively. Comparison between the curves facilitates 
the identification of the incident channel responsible for the 
Feshbach resonances. Magnetic field resolution is lO/iT. 



E. (3,3) + (4,2), (3,3) + (4,3), and (3,3) + (4,4) 
inelastic collisions 

Multiple Feshbach resonances due to (3, 3) + (4, 2) and 
(3, 3) + (4, 3) scattering are observed in samples with 80% 
population in |3,3), and 15% in |4, 3) or |4, 2). The re- 
maining atoms are predominantly in the |3, 2) hyperfine 
state. We have verified that (3, 3) + (3, 2), (4, 3) + (3, 2) 
and (4, 2) + (3, 2) processes do not contribute to the loss. 

We prepare the sample by first polarizing 95% of the 
atoms in the |3,3) state at _B = 14/iT, and then apply a 
microwave pulse for a few ms, which selectively transfers 
approximately 15% of the population into |4,3) or |4, 2). 
Due to the almost identical energy splitting between the 
3,3) ^ |4, 2) and |3,2) ^ |4, 3) microwave transitions 
at low field, we have a small population of |4, 3) in the 
experiment aimed at finding resonances in a (3, 3) + (4, 2) 
collision. Pollution by |3, 3} + |4, 3) resonances occurs, as 
was previously reported in Ref. p?7| . 

We also observe a weak pollution by (3, 3) + (4, 2) res- 
onances in the (3, 3) -I- (4, 3) spectrum. This cannot be 
explained by microwave transitions, since |3, 3) ^ |4, 3) 
and |3, 2) ^ |4, 2) do not have the same frequency. A pos- 
sible process that creates atoms in the |4, 2) state is the 
inelastic collision process (3, 2) -|- (4, 3) (3, 3) -I- (4, 2). 
At low fields the hyperfine and Zeeman energy difference 
between the final and initial state is 6E = fcB0.52/iK 
x{B/mT). For our atomic temperatures and magnetic 
fields, this endothermal spin-changing collisions can cre- 
ate atoms in the |4, 2) state. 

The four resonances in the (3, 3) + (4, 2) channel at 
6.17 mT, 8.38 mT, 11.0 mT and 11.2 mT, and all four 
resonances observed in the (3, 3) -I- (4, 3) channel at 12.9 
mT, 17.3 mT, 22.7 mT and 23.0 mT are paired with 
identical ratio of the field values 2.07(2). This observa- 
tion is confirmed by the theoretical identification of the 
paired resonances as being due to molecular states with 
identical binding energy at zero magnetic field. 

In a separate experiment we have found no inelastic 
resonances in (3, 3) -|- (4, 4) collisions for magnetic field 
strengths up to 23.5 mT, in agreement with theory. 



IV. ELASTIC FESHBACH SPECTROSCOPY 



D. (3,3) + (3,3) and (3,3) + (3,2) inelastic collisions 

The (3, 3) -I- (3, 3) channel is the lowest hyperfine scat- 
tering channel, and therefore has no binary exothermic 
collision processes. Collisional inelastic loss is then only 
due to three-body recombination, a process that falls out- 
side the scope of this paper. A quantitative study of the 
recombination loss in the |3, 3) state is given in Ref. [28|| . 
In collisions between |3,3) and |3,2) atoms we have ob- 
served no resonances below 5=23.5 mT. Theory predicts 
weak resonances that are beyond the sensitivity of our 
current experiment. 



The resonant change of the scattering amplitude in the 
incident channel also results in a modification of the elas- 
tic cross section. We refer to this process as an elastic Fes- 
hbach resonance, see Fig. ^ Beside the direct measure- 
ment of the thermalization rate between the axial and ra- 
dial motion [isL |2^ , we have developed a more sensitive 
measurement technique for elastic Feshbach resonances, 
the measurement of the evaporation rate in a shallow 
trap 01 . This method converts the temperature evo- 
lution measurement into an atom number measurement, 
which provides better sensitivity and signal-to-noise ra- 
tio. 

We report elastic collision properties in the pure 
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(3,3) + (3,3) channel and in the mixed (3,3) + (3,2) 
channel, where the measurements are not complicated 
by inelastic processes. The elastic spectrum of the 
(3, 3) + (3, 3) channel has previously been reported in 
Refs. [l7l [l8l|. Although the elastic cross section can in 
principle be both enhanced and decreased by the Fesh- 
bach resonance, we observe only dips in the elastic cross 
section at T=5 fiK. This is due to the large background 
scattering length, resulting in an elastic cross section that 
is mostly unitarity hmited to ct = Stt/A:^ (ct = 47r/A:^) 
when the colliding cesium atoms are (are not) in the same 
state. A vanishing elastic collision cross section at 1.7 
mT resulting from a strong s— wave Feshbach resonance 
in the (3, —3) + (3, —3) channel was first observed as a 
reduction in the cross axis thermalization fl^ . Three ad- 
ditional c?— wave Feshbach resonances were observed by 
monitoring the evaporative loss rate ^ vanishing 

elastic cross section in the (3, 3) + (3, 2) channel is ob- 
served with the same evaporative loss method, where we 
measure the remaining atom number in |3,2). Positions 
of the elastic collision rate minima are given in Table |I| 
The elastic collision measurements based on evapora- 
tive loss are limited to resonances whose width in energy 
is comparable to or exceeds the thermal energy of the 
sample kgT. For narrow resonances the resonant varia- 
tion of the elastic collision rate over a small interval of 
collision energies does not lead to observable variations 
in the evaporation rate. For the detection of such weak 
resonances we have to resort to radiative Feshbach spec- 
troscopy, discussed in the next section. 



V. RADIATIVE FESHBACH SPECTROSCOPY 

Narrow resonances with a linewidth that is small com- 
pared to the thermal energy ksT do not significantly af- 
fect the sample's thermalization rate, and therefore can- 
not be probed by elastic Feshbach spectroscopy. How- 
ever, it is possible to directly measure the enhanced 
(quasi-)bound state population on a Feshbach resonance 
using a far-detuned laser beam to selectively dissoci- 
ate the molecules, while leaving the atoms unperturbed. 
This technique is called radiative Feshbach spectroscopy. 

In this work, we apply radiative Feshbach spectroscopy 
to probe Feshbach resonances in the (3, 3) + (3, 3) and 
(3,3) -|- (3,2) scattering channels using a probe beam 
typically detuned 20 GHz to 4 THz to the blue of the 
cesium D2 transition at 852.3 nm. The probe beam is 
provided by a titanium-sapphire laser, and uniformly il- 
luminates the atom sample with a stabilized intensity up 
to 50W/cm^. The intensity and detuning of the laser 
are adjusted such that the single-atom excitation is suf- 
ficiently weak, and negligible atom loss is observed when 
the magnetic field is tuned off the Feshbach resonances, 
while maximizing the loss on the Feshbach resonance. 

To measure the radiative collision loss, we first pre- 
pare atom samples either fully polarized in the |3, 3) to 
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FIG. 6: A detailed radiative loss spectrum with a field reso- 
lution of 0.5/iT. The probe beam at wavelength A — 844 nm 
has an intensity of 50W/cm^. 85% (15%) of the atoms are in 
the |3, 3) (|3, 2)) state. The temperature and mean density of 
the sample are 3.5/^K and 1 x 10^'^ cm~^, respectively. "S" 
indicates a shape resonance; the stars indicate Feshbach reso- 
nances in the (3, 3) -I- (3, 2) channel. Magnetic field resolution 
near the resonances is 5/iT. 



study (3, 3) + (3, 3) colhsions, or 85% in |3, 3) and 15% in 
|3, 2} to study (3, 3) -t- (3, 2) collisions. After the prepara- 
tion and field ramping, we typically illuminate the atoms 
with the probe beam for t=100 ms to 300 ms, and mea- 
sure the trap loss. To determine the radiative loss rate as 
a function of magnetic field, we perform two consecutive 
atom number measurements at every magnetic field value 
B = Bi to discriminate non-radiative loss: Ni is the atom 
number at B = Bi without the probe beam, and N2 is 
obtained at B = Bi with the probe beam. The radiative 
loss rate is then given by 7 = {1 — N2/Ni)/t, where the il- 
lumination time t is chosen to keep the maximum atomic 
loss below 30%. To increase detection sensitivity on weak 
resonances, we lengthen t to > 500ms, and average each 
data point up to 5 times. A detailed spectrum thus ob- 
tained is shown in Fig.|Bl The details of the radiative loss 
line shape are discussed in Ref. 0|; atom-molecule dy- 
namics and the sensitivity of the our radiative Feshbach 
spectroscopy are studied in Ref. . Positions of the res- 
onances are tabulated in Table IIIII One of the Feshbach 
resonances observed by radiative Feshbach spectroscopy 
was recently used to create ultracold cesium molecules 



i 



VI. THEORY 

The structure of the Hamiltonian of two interacting 
^S* ground state alkali- metal atoms is well known. It 
contains the atomic kinetic energy operator, an atomic 
Hamiltonian for each atom, two Born-Oppcnheimer po- 
tentials with symmetry and '^Sj, the nuclear rota- 
tion operator h^P / {2^R^), and weaker relativistic spin- 
spin dipole and second order spin-orbit interactions. Here 
I is the nuclear mechanical angular momentum and /i is 
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the reduced mass of the molecule. 

The atomic Hamiltonian contains a Fermi contact term 
and the Zeeman interaction when an external magnetic 
field, B, directed along the z axis is present during the 
collision. The eigenstates of the B — atomic Hamilto- 
nian are jFo,, rnf^a), where a = 1 or 2 for atom 1 or 2, 
Fa = Sa + ia and mpa is the projection of Fa along the 
z axis. Here Sa and ia are the atomic electron and nu- 
clear spin, respectively. For B > 0, states with the same 
mpa mix, Fa is no longer a good quantum number, and 
the zero field mpa degeneracy is lifted. For convenience, 
atomic eigenstates in a magnetic field will be labeled by 
iFa^mpa) since for fields used in this paper the Zeeman 
interaction is small compared to the hyperfine interac- 
tion and Fa is an approximately good quantum number. 
We use hyperfine constants and magnetic moments from 
Ref. [iOl. The nuclear spin of cesium is 7/2. 

The selection rules for the Born-Oppenheimer Hamil- 
tonian conserve / and f = Fi + F2. Consequently, bound 
states and scattering amplitudes can be labeled by fl for 
B = and rrifl for B > 0. Here, mf — nipi + mp2. 
The two relativistic interactions weakly mix states with 
different I and /. Global symmetries ensure that the 
molecular Hamiltonian, including the relativistic inter- 
actions, conserves parity and total angular momentum 
F = f + I and its projection M = mf + mi for i? = 
or only M for nonzero fields. Only even or odd partial 
waves are coupled. For B ^ there are at most 72 cou- 
pled channels while for non-zero field there are infinitely 
many coupled channels. In practice the number of chan- 
nels is restricted by using knowledge about the relative 
strengths of the individual terms in the Hamiltonian. 

The relativistic interactions, even though weak, are 
crucial in understanding the presence of 5- wave Feshbach 
resonances in our gas of ultra-cold Cs atoms. The tem- 
perature T « 5 /iK is small compared to the 200 /iK and 
1 mK barrier height of the d- and g-partial waves, respec- 
tively, and only incoming s-wave collisions contribute to 
the experimental signal. Consequently, coupling from s- 
to d- and 5- waves and back to s-waves is necessary. Small 
corrections to the location of the Feshbach resonances 
due to relativistic interactions can be observed. 



VII. NUMERICAL APPROACHES 

The scattering properties and bound-state energies of 
the ground state Hamiltonian are obtained with two sep- 
arate numerical approaches. For scattering the wavefunc- 
tion at energy E is found using a Gordon propagator |3l|. 
From the scattering wavefunctions elastic and inelastic 
rate coefficients as a function of magnetic field can be 
evaluated. Feshbach resonances appear as sharp peaks 
or dips in the magnetic field dependence of the rate co- 
efficients. For comparison with the experiments the rate 
coefficients need to be thermally averaged. 

Obtaining discrete bound states with the Gordon 
method is cumbersome since eigenenergies are not a pri- 



ori known and multi-channel scattering wavefunctions 
need to be calculated at a large number of energies E. 
Consequently, a discrete variable representation [22| for 
the radial kinetic energy operator is used to find the 
bound states. In this approach the eigenvalues of a lin- 
ear system of size given by the number of radial colloca- 
tion points times the number of coupled channels need 
to be calculated. This can be done with standard lin- 
ear algebra packages. However, resource limitations tend 
to restrict the number of coupled channels that can be 
conveniently handled. For the heavy cesium dimer a re- 
alistic maximum number of channels lies between 10 to 
15, although 20 channels can still be treated. An inter- 
esting alternative for finding bound state energies, which 
does not require the storage of the linear system, can be 
based on the multi-channel quantum defect theory [33j |. 
It takes advantage of analytic properties of wavefunctions 
as a function of energy E in order to limit the number of 
wavefunction evaluations. 

Bound states are calculated over a range of magnetic 
field values. Feshbach resonances occur when a bound 
state crosses a collisional threshold. Examples of reso- 
nances can be found in Refs. 0) Ull' In the absence 
of the two relativistic interactions calculations label in- 
dividual bound state by to/, /, and M. In addition, the 
field dependence of a level can be traced to a bound state 
at zero field. At zero field / labels the bound states. For 
magnetic fields used in our experiment / remains approx- 
imately good. Coupling of states with different mf but 
the same I are sometimes needed to fully assign the Fesh- 
bach resonances. For this paper bound state calculations 
are used to assign quantum numbers to the resonances. 

The Feshbach resonances are experimentally observed 
in either elastic cross sections, inelastic rate coefficients, 
or radiative collision rates. The former two measure- 
ments can be modeled from first-principle scattering cal- 
culations of cross sections and rate coefficients. In prin- 
ciple the rate coefficients need to be thermally averaged. 
However, a proper thermalization was impractical and 
for a comparison between theory and experiment only a 
single collision energy given by the mean collision en- 
ergy of a gas at temperature T was used. A combi- 
nation of narrow Feshbach resonances and the need to 
study the effect of variations in the shape of the two 
Born-Oppenheimer potentials on the resonance locations 
would lead to an untenable number of scattering calcula- 
tions. A one-standard-deviation uncertainty of 0.02 mT 
to 0.05 mT in the calculated magnetic-field location of 
the Feshbach resonances observed in elastic or inelastic 
rates results from the use of a single collision energy. 

The radiative collision data require modeling of the 
rate coefficient for the absorptions of a photon by a pair 
of ultra-cold Cs atoms ^36). However, we are interested 
in reproducing the location of Feshbach resonances, and 
not in the absolute absorption rates. Consequently, from 
a theoretical perspective it is sufficient to locate the reso- 
nances in the elastic cross section in the absence of light. 
The I = 4 ( (7- wave) resonances observed by radiative Fes- 
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hbach spectroscopy could not be observed by direct mea- 
surement of the elastic cross section. 

The energy widths oi I — A (g-wave) Feshbach reso- 
nances are significantly smaller than the average collision 
energy. As a consequence, the location of the maximum 
of the photoassociation lineshape as a function of mag- 
netic field depends on the zero-collision-energy resonance 
field location Bq, the temperature, and the magnetic mo- 
ment /Ltres of the embedded bound state. In practice, how- 
ever, we located the resonances from scattering calcula- 
tions at E/kg — 5.3 //K and used the magnetic moment 
of the (7- wave Feshbach resonance obtained from bound- 
state calculations to extrapolate to zero collision energy. 



VIII. RESULTS 

Tables ^1 ^nd IIIII give the magnetic field locations 
and assignments of the observed Feshbach resonances. 
Numbers in parentheses indicate the one standard de- 
viation uncertainty. The locations of the resonances are 
obtained from elastic and inelastic cross-section measure- 
ments or radiative spectroscopy. The theoretical reso- 
nance locations are obtained from coupled-channel scat- 
tering calculations at a E/kB=5.3 /xK collision energy 
and Born-Oppenheimer potentials with dispersion coeffi- 
cients 6*6=6890 EhttQ and Cs =954600 Eha^, and scatter- 
ing lengths as = 280.3 uq and ax = 2405 gq |15ij. Here 
ks is the Boltzmann constant, 1 Eh — 4.35974 aJ is a 
Hartree, and Iuq = 0.0529177 nm is a Bohr. kg^.S fiK 
is the average collision energy for a Cs gas at T=3.5 ^K. 
For resonances observed in the elastic (inelastic) rates the 
minimum (maximum) of the line is quoted. The only ex- 
ception is the resonance at 7.66 mT, where the position of 
the minimum inelastic rate is given. For the resonances 
observed by the radiative spectroscopy the maximum of 
the loss rate and the maximum of the theoretical inelastic 
collision rate in the absence of laser light are presented. 

It should be noted that the magnetic field values for 
the spectral features presented here should not be con- 
fused with the location of the molecular Feshbach state 
introduced in Fig. ^ Typically, the molecular state can 
be defined anywhere within the width of the resonance. 

The theoretical uncertainties are a consequence of our 
limited ability to model the experiments, and are ob- 
tained by combining the uncertainties due to the spread 
in collision energies of a thermal Cs gas with the mag- 
netic field dependence of the Feshbach resonance. Any 
discrepancy between theory and experiment in the tables 
that lies outside the error bars indicate deficiencies in the 
shape of the two Born-Oppenheimer potentials. 

The assignment of the resonances is obtained by com- 
bining information from scattering and bound state cal- 
culations. The initial collision partners that lead to a res- 
onance are determined from theoretical scattering calcu- 
lations, and experimentally by varying the relative popu- 
lation of hyperfine states in the Cs gas and comparing the 
relative strength of the resonances. The incoming par- 



TABLE I: Location and assignment of the observed Feslibacli 
resonances. The first three columns define and give the results 
of the experiment, where the first column denotes the initial 
collision state, the second column indicates elastic (el.), in- 
elastic (inel.), or radiative (rad.) measurements, and the third 
column is the experimental resonance location Bexp- The next 
four columns describe the quantum labels of the resonances. 
These are the partial wave of the initial collision state, the 
partial wave of the Feshbach state, the molecular spin /, and 
its projection m obtained from multichannel bound state cal- 
culations. The last column gives the theoretical resonance 
location Bth obtained from a multichannel scattering calcula- 
tion at a collision energy of E/kB=5.3 fiK. The resonance at 
Bex ~ 20.66 mT was not observed in the theoretical scatter- 
ing calculation due to its extreme narrowness but nevertheless 
could be assigned from bound state calculations. Several Fes- 
hbach resonances have ambiguous assignments (see Text). 
Experiment Theory 

State Method i?..(mT) ' Assigniiien^ i?,,(mT) 

Wave I f ruf 



(3,3) 


+ 


(3,3) 


el. 


^1.706(3) 


s 


s 


6 


6 


1.70(2) 


(3,3) 


+ 


(3,3) 


el. 


^4.802(3) 


s 


d 


4 


4 


4.79(2) 


(3,3) 


+ 


(3,2) 


el. 


^5.69(2) 


s 


d 


4 


4 


5.70(2) 


(4,-4) 


+ 


(4,-4) 


inel. 


10.590(3) 


s 


d 


8 


-6 


10.58(2) 


(4,-4) 


+ 


(4,-4) 


inel. 


20.503(3) 


s 


d 


8 


-7 


20.49(2) 


(4,-4) 


+ 


(4,-3) 


inel. 


20.66(1) 


s 


d 


8 


-6 




(3,-3) 


+ 


(3,-3) 


inel. 


3.005(5) 


s 


d 


6 


Ri-6 


2.99(2) 


(3,-3) 


+ 


(3,-3) 


inel. 


3.305(5) 


s 


d 


6 


-4 


3.28(2) 


(3,-3) 


+ 


(3,-3) 


inel. 


8.69(2) 


s 


d 


8 


-8 


'^8.80(2) 


(3,-3) 


+ 


(3,-3) 


inel. 


10.11(2) 


s 


d 


8 


-7 


10.15(2) 


(3,-3) 


+ 


(3,-3) 


inel. 


10.88(2) 


s 


d 


8 


-6 


10.90(2) 


(3,-3) 


+ 


(3,-3) 


inel. 


11.81(2) 


s 


d 


8 


-5 


11.85(2) 


(3,-3) 


+ 


(3,-3) 


inel. 


13.31(2) 


s 


d 


8 


-4 


13.35(2) 


(3,3) 


+ 


(4,2) 


inel. 


6.17(2) 


s 


d 


7 


5 


6.21(2) 


(3,3) 


+ 


(4,2) 


inel. 


■^7.66(2) 


s 


s 


6 


5 


7.53(2) 


(3,3) 


+ 


(4,2) 


inel. 


8.05(2) 


p 


P 


6 


5 


8.08(2) 


(3,3) 


+ 


(4,2) 


inel. 


8.38(2) 


s 


d 


7 


6 


8.43(2) 


(3,3) 


+ 


(4,2) 


inel. 


11.00(3) 


s 


d 


7 


5,7 


11.02(2) 


(3,3) 


+ 


(4,2) 


inel. 


11.20(3) 


s 


d 


7 


5,7 


11.20(2) 


(3,3) 


+ 


(4,2) 


inel. 


16.22(4) 


s 


s 


5 


5 


16.23(2) 


(3,3) 


+ 


(4,2) 


inel. 


18.29(5) 


p 


P 


5 


5 


18.45(2) 


(3,3) 


+ 


(4,3) 


inel. 


12.90(3) 


s 


d 


7 


6 


12.96(2) 


(3,3) 


+ 


(4,3) 


inel. 


17.30(4) 


p 


P 


6 


6 


17.45(2) 


(3,3) 


+ 


(4,3) 


inel. 


22.73(5) 


s 


d 


7 


5,7 


22.75(2) 


(3,3) 


+ 


(4,3) 


inel. 


23.05(5) 


s 


d 


7 


5,7 


23.13(2) 



" Temperature-dependent resonance occurred at Re{a) = 0. 
^ Minimum evaporation rate. 
* Minimum inelastic loss. 
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TABLE II: Feshbach resonances between two F = 3 atoms ob- 
served in measurements of inelastic rate coefficients. Columns 
defined as in Table One bound state could not be assigned 
while three resonances have not been observed experimentally. 

Experiment Theory 





State 


Method 


Be.(mT) 


Inc. 
Wave 


Assij 

I f 




Bth(mT) 


(3,-3)+(3,-3) 


inel 


2.18(2) 










* 


(3, 


-3)+(3,-2) 


inel. 


3.57(2) 


s 


d 


6 


-5 


3.55(5) 


(3, 


-3)+(3,-2) 




1 n 50^ 1 


P 


f 


7 


_7 


in 5n('5') 


(3, 


-3)+(3,-2) 


inel 


11.04(2) 


P 


f 


7 


_7 


11.03(5) 


(3, 


-3)+(3,-2) 


inel 




P 


f 


7 


7 




(3, 


-3)+(3,-2) 


inel. 


12.01(2) 


s 


d 


8 


-7 


12.03(5) 


(3, 


-3)+(3,-2) 


inel. 


13.01(2) 


s 


d 


8 


-6 


13.05(5) 


(3, 


-3)+(3,-2) 


mel. 




s 


d 


8 


-o 


1 A 


(3, 


-3)+(3,-2) 


inel. 




s 


d 


8 


A 


1 7 

1 / .uo^^o ) 


(3, 


-3)+(3,-2) 


inel. 




s 


d 


8 


-3 


'20.88(5) 


(3, 


-2)+(3,-2) 


inel. 


^14.82(2) 


s 


d 


8 


-7 


14.80(5) 


(3, 


-2)+(3,-2) 


inel. 


^16.58(2) 


s 


d 


8 


-6 


16.50(5) 


(3, 


-2) + (3,-2) 


inel. 


^19.25(2) 


s 


d 


8 


-5 


19.25(5) 


(3, 


-2) + (3,-2) 


inel. 




s 


d 


8 


-4 


■23.75(5) 


(3, 


-3+3,-1) 


inel. 


12.92(2) 


s 


d 


8 


-8 


12.90(5) 


(3, 


-3+3,-1) 


inel. 




s 


d 


8 


-4 


■21.15(5) 



* not predicted by calculation. 
' not observed experimentally. 

^ has equal contribution from (3,-3) + (3,-l) collision. 



tial -wave is obtained from scattering calculations. The 
assignment of quantum numbers /, /, and is made 
on the basis of bound-state calculations. One resonance 
could not be assigned. 

Levels "with the same /,/,m^ symmetry but different 
M = mi + m f are degenerate except for small splittings 
from the second-order spin-orbit and spin-spin interac- 
tions. For even partial waves this is not an issue as losses 
from the s-"wave entrance channel at collision energies of 
the order of a fe^w microkelvin are much larger than those 
from d-wave entrance channels and thus only resonances 
with M equal to the sum of the magnetic quantum num- 
ber of the initial hyperfine states, M = rufa + nifb, are 
observed. For collisions between atoms in unlike hyper- 
fine states contributions from p-wave collisions cannot be 
ignored and nearly degenerate Feshbach resonances from 
three M values appear. 

A good example of the complexity for odd I resonances 
can be found between 4.0 mT and 6.0 mT in the radiative 
Feshbach spectrum from a (3, 3) + (3, 2) collision. In this 
region six nearly degenerate resonances are labeled / ~ /, 
f = 5, ruf — 2 or I = p, f = 5. The m/ labels of the p- 
wave resonance could not be assigned. It turns out that 
this resonance has a small magnetic moment and at _B w 



TABLE III: Feshbach resonances and a single shape resonance 
observed with radiative spectroscopy. Columns defined as in 
Table Several Feshbach resonances have ambiguous assign- 
ments (see Text). 



Experiment Theory 



State 


Method 


Be.(mT) 


Inc. 
Wave 


Assignment 
/ / nif 


Bth(mT) 


(3,3)+(3,3) 


rad. 


1.102(3) 


s 


g 4 


2 


1.12(2) 


(3,3)+(3,3) 


rad. 


1.437(3) 


s 


g 4 


3 


1.46(2) 


(3,3)+(3,3) 


rad. 


1.506(3) 


s 


g 6 


5 


1.51(2) 


(3,3)+(3,3) 


rad. 


1.83(1) 


s 


d 




n.86(2) 


(3,3)+(3,3) 


rad. 


1.990(3) 


s 


g 4 


4 


2.01(2) 


(3,3)+(3,3) 


rad. 


4.797(3) 


s 


d 4 


4 


4.77(2) 


(3,3)+(3,3) 


rad. 


5.350(3) 


s 


g 2 


2 


5.43(2) 


(3,3)+(3,3) 


rad. 


11.278(3) 


s 


d 6 


4 


11.32(2) 


(3,3)+(3,3) 


rad. 


13.106(3) 


s 


d 4 


4 


13.19(2) 


(3,3)+(3,2) 


rad. 


0.78(1) 


s 


g 6 


3 


0.83(2) 


(3,3)+(3,2) 


rad. 


1.13(1) 


s 


g 4 


1 


1.17(2) 


(3,3)+(3,2) 


rad. 


1.47(1) 


s 


g 4 2 or 6 4 


1.54(2) 


(3,3)+(3,2) 


rad. 


1.66(1) 


p 


f 3 


2 


1.64(2) 


(3,3)+(3,2) 


rad. 


2.09(1) 


s 


g 4 


3 


2.16(2) 


(3,3)+(3,2) 


rad. 


2.21(1) 


p 


f 1 


1 or 3 3 


2.18(2) 


(3,3)+(3,2) 


rad. 


2.36(1) 


p 


f 1 


1 or 3 3 


2.33(2) 


(3,3)+(3,2) 


rad. 


3.60(1) 


s 


g 4 


4 


3.70(2) 


(3,3)+(3,2) 


rad. 


3.81(1) 


P 


f 5 


1 


3.78(2) 


(3,3)+(3,2) 


rad. 


4.68(1) 


p 


f 5 2 


or p /=5 


4.70(2) 


(3,3)+(3,2) 


rad. 


4.93(1) 


p 


f 5 2 


or p /=5 


4.89(2) 


(3,3)+(3,2) 


rad. 


4.99(1) 


p 


f 5 2 


or p /=5 


4.99(2) 


(3,3)+(3,2) 


rad. 


5.70(1) 


s 


d 4 


4 


5.70(2) 


(3,3)+(3,2) 


rad. 


5.77(1) 


p 


P 5 


3,4,5 


5.78(2) 


(3,3)+(3,2) 


rad. 


5.87(1) 


p 


P 5 


3,4,5 


5.86(2) 


(3,3)+(3,2) 


rad. 


5.97(1) 


p 


P 5 


3,4,5 


5.98(2) 



d-wave shape resonance 



5 mT the different m / components are nearly degenerate. 
Mixing between different m f components for the same M 
due to second-order spin-orbit and spin-spin interactions 
leads to shifts that are comparable to the spacings due 
to the Zeeman interaction. The / = / resonances, which 
have a larger magnetic moment, accidentally reside in the 
same magnetic field region. 

We have assigned more than one set of labels to the 
resonances between 4.0 mT and 6.0 mT, because m/ and 
to a lesser degree I and / are not good quantum labels. 
There are eight different I = p, f = 5 and I = f , f — 5, 
nif — 2 resonances between 4.0 mT and 6.0 mT. Six of 
these eight resonances are due to the p-wave symmetry. 
This can easily be checked by noting that only M—A, 5, 
and 6 can lead to these odd I resonances. Not all eight 
resonances have been seen which might be due to the 
fact some are not resolved or unobserved by radiative 
spectroscopy. We did not perform quantitative bound 



/ 4d 

-0 15 1 1 1 ' 1 1 1 1 1 1 1 1 1 1 1 

5 10 15 

B (mT) 

FIG. 7: Total angular momentum projection M=6, s, d, and 
jr-wave bound state energies as a function of magnetic field. 
The zero of energy corresponds to the (F = 3, mp — 3) + {F = 
3, m_F = 3) dissociation limit. Dotted, dashed, and solid lines 
correspond to l=s,d, and g states respectively. Furthermore, 
each curve is labeled by the quantum numbers fl, i.e. fl — 6s. 
The molecular spin / is a zero field quantum number. For 
the fl = 6d states the magnetic quantum number m/ is also 
indicated. The filled circles represent the observed threshold 
resonances. 



state calculations as too many channels must be coupled 
together. The theoretical field locations listed in the last 
column of the tables have been obtained from scattering 
calculations where all states are included at the cost of 
losing the ability to assign quantum labels. 

Fig. H shows even I M=6 bound states below the low- 
est molecular hyperfine state (3,3) -I- (3,3) as a function 
of magnetic field. Each bound state is labeled with I, f, 
and nif. A resonance occurs when a bound state crosses 
zero energy. The frequency range shown in the figure is 
sufficient for the assignment of all i? <15 mT s, d, and g- 
wave Feshbach resonances in the collision between |3,3) 
Cs atoms. The filled circles mark the observed thresh- 
old resonances in a Cs gas at T=3.5 /xK {ksT/h ~ 0.1 
MHz) . Agreement between theory and experiment is suf- 
ficiently good that assignments can be made even though 
discrepancies exist. These discrepancies are caused by 
the (slightly) incorrect shape of the Born-Oppenheimer 
potentials and the approximations in the bound state cal- 
culations. For fields smaller than 1 mT theory predicts 
the existence of additional resonances. 

The number of coupled channels for A/=6 and l—s^ 
d, and g is 74. However, as discussed in the previous 
section, for nonzero applied magnetic field the coupling 
between different partial waves and m / 's is due to weak 
second-order spin-orbit and spin-spin interactions. Con- 
sequently, for most resonances in Fig. [7| assignment is 
unambiguous using independent bound-state calculations 
that only include states of a given / and toj. In fact, 
the curves in the figure have been obtained in this way. 
However, it should be realized that crossings between 





1.4 1.45 1.5 1.55 



B (mT) 

FIG. 8: Expanded view of the M=6, s, d, and g-wave bound 
states shown in Fig. Q An avoided crossing between g-wa,ve 
mf—3 and 5 bound state occurs around B= 1.4 mT. 



bound states shown in the figure are actually avoided 
when second-order spin-orbit and spin-spin interactions 
are included. At zero magnetic field coupling between 
channels with different molecular spin / is also small. 
The assignment of / is obtained by retracing a bound 
state to zero magnetic field and noting that < /. 

A close look at Fig. |7| shows that the lines can roughly 
be divided into those that are noticeably curved and 
those that appear straight. A good example of curved 
lines are the two fl = 4d bound states, while the 6d and 
Ag bound states are good examples of bound states that 
have a linear magnetic field dependence. The curved lines 
are due to broad avoided crossings that appear because 
the same fl labeled states readily mix when a magnetic 
field is applied. Mixing is due to the interplay of hyper- 
fine, Born-Oppenheimer, and Zeeman interactions and is 
significantly larger than in avoided crossings mediated by 
second-order spin-orbit and spin-spin interactions. 

The most weakly-bound fl — 6s state is bound by 
about 65 MHz at zero magnetic field, rises rapidly until 
it turns over near 2 mT, and then continues just be- 
low the dissociation limit. This avoided crossing is also 
shown in Fig. |S1 The bound state does not run parallel 
to the dissociation limit. It becomes a Feshbach reso- 
nance near B—50 mT. The behavior of this bound state 
has direct consequences for the s-wave scattering length 
of two IfaiTia) — |3,3) atoms. Below 1.7 mT the scat- 
tering length is negative and above this field value it is 
positive. This zero of the scattering length has been ob- 
served in Ref. jl3 and has been used t o op timize the 
Born-Oppenheimer potentials in Refs [T^ [13 , as well as 
in this paper. It is interesting to realize that, as discussed 
in Ref. [l^, for B < 1.7 mT d-wave channels affect the 
elastic scattering and must be included in order to obtain 
an accurate scattering length. 

Some of the M=6 and ^=s, d and g Feshbach reso- 
nances below B=3 mT could at first not be assigned from 
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FIG. 9; Position of two M = 6 g-wsye resonances as a func- 
tion of Cg or equivalently Ce for a zero-energy collision of two 
Cs atoms in the lowest hyperfine state. The Ce and Cg are 
expressed in units of E^aQ and E^aQ, respectively. The dot- 
ted lines correspond to the peak of the experimental radiative 
Feshbach spectroscopy signal. 




5 10 15 

magnetic field 6 [mT] 



FIG. 10: Scattering length in the {F = 3, = 3) -|- (F = 
3, mF =3) scattering channel. Resonances resulted from the 
M=6, s, d, and g-wave bound states are assigned with quan- 
tum number (Z, f,mj), where I is the orbital angular momen- 
tum, / is the total internal angular momentum and m/ is 
the magnetic quantum number. Calculation are done with a 
magnetic field grid size of 50/iT for off-resonance regions and 
lOOnT near the narrow resonances. Inset shows the detail 
resonance structure near 17G 



calculations using states with the same I and nif. Res- 
onances of different Im j symmetry lie in the same mag- 
netic field region. Fig. |S1 shows a blow up of the mT 
to 3 mT magnetic field range. The /=4 m/=3 and /=6 
m/=5 g-wave bound states cross just below the dissocia- 
tion hmit and weak couplings might shift the correspond- 
ing Feshbach resonances. The inset shows the avoided 
crossing between these /=4 mf—3 and /=6 mf—5 g- 
wave bound state when the weak couphng between the 
two bound states is included. From the figure it is clear 
that the avoided crossing has little influence on the lo- 
cation of the Feshbach resonances and an fm / label for 
each resonance can still be assigned. 

The Born-Oppenheimer potentials that have been used 
for the theoretical resonance locations quoted in TablesQl 
HH andHniand Figs. Hand El are based on the fit in our 
previous work [15|, where the major uncertainties in the 
calculation of resonance positions arise from the poorly 
constrained Cs coefficient. For this Paper we have im- 
proved the Born-Oppenheimer potentials by optimizing 
the Cs dispersion coefficient in addition to the Cg coef- 
ficient, the strength of the second-order spin-orbit inter- 
action Sc, and the singlet and triplet scattering lengths 
[Tsj l . For a given Cs the potentials are optimized to fit the 
minima in the elastic scattering rate of the (3, 3) -I- (3, 3) 
scattering at 1.7064 mT and 4.8017 mT, the / = 6 d- 
wave resonance in (3, —3) -I- (3, —3) scattering at 3.0051 
mT, and the resonance in (4, —4) -I- (4, —4) scattering at 
20.5029 mT. Over a 10% range of the Cs coefficient near 
900000 EhttQ a linear relationship between Cg, Sc, as, 
and gt exists. 

An improved cesium-dimcr Hamiltonian is created by 
fitting to selected (3, 3) -f (3, 3) g-wa,ve resonances. Fig. El 
shows the M=6, /=6 mf=5 and /=4 m/=3, g-wave res- 
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magnetic field B [ml] 

FIG. 11: Collision rate constants in the {F — 3,mF = 
—3) + {F — 3, mp ~ —3) scattering channel at E/ks ~ InK. 
Inelastic collision rates (a) and elastic collision rates (b) are 
calculated with a magnetic field grid size of lO/xT for off- 
resonance regions and l/xT near the resonances. Resonances 
from the M=-6, s, d, and g-wave bound states are included 
and assigned with quantum number {l,f,mf), where the no- 
tation is the same as in Fig. IIUI 



onance position as a function of Cs or equivalently Cq. 
The figure shows theoretical zero-collision-energy reso- 
nance locations derived from E/kB=5.3 /iK calculations 
and the magnetic moments of the resonances. The peak 
radiative detection signal as a function of magnetic field 
is due to collisions at zero energy. The magnetic moment 
of the resonances is 170 /xK/mT for the /=6, to/=5 and 
550 AtK/mT for the /=4, m/=3 state. 

The location of the Feshbach resonance found from a 
bound state calculation, including only g-wave mf=3 and 
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TABLE IV: Properties of the singlet X^E+ and triplet a^Sj 
Born-Oppenheimer potentials and the second-order spin-orbit 
interaction that give the best fit to all data on collisions be- 
tween ultra-cold Cs atoms. The Ce and Cg are expressed in 
units of Ehag and EnaQ, respectively. The singlet, as, and 
triplet, At, scattering lengths are in units of oq. Sc is dimen- 
sionless. One-standard deviation uncertainties are given. 





value 


uncertainty (%) 




6860 


0.36 


Cs 


860 000 


8.7 


as 


280.37 


0.02 




2440 


1.0 


Sc 


2.6 


19 



5 channels, and a scattering calculation at zero collision 
energy, which includes all s-, d-, and g-wave channels, 
disagree by about 0.02 mT. This discrepancy is likely 
due to the limited number of channels in the bound state 
calculations. The magnetic moments, however, are not 
expected to be significantly modified. 

Table summarizes our best fit. Based on the colli- 
sion parameters, cold collision properties of cesium atoms 
can be readily calculated in various scattering channels. 
In particular, (3, 3)-|-(3, 3) and (3, — 3)-|-(3, —3) scattering 
length is —2510 ao at zero magnetic field. In the presence 
of the magnetic field. Fig. [TUl shows the scattering length 
of the (3,3) -I- (3,3) channel as a function of magnetic 
field; Fig. ^] shows the collision rate constants in the 
(3, — 3)-|-(3, —3) channel at the collision energy oi E — 
InK. All s, d, and g wave channels are included in the 
calculation. These two states are particularly interesting 
in the experiments of cesium Bose-Einstein condensation 
[37| . The resonances are identified and labelled by the 
quantum numbers of the associated molecular states. 

IX. CONCLUSION 

We have measured > 60 magnetic-field-induced Fesh- 
bach resonances in the collision of ultracold ground state 
cesium atoms. Of all the alkali-metal species cesium is 
shown to have the richest resonance structure. 

The field resonances have been observed in elastic col- 



lision rates, evaporation loss rates, collision relaxation 
rates, as well as in radiative collision resonance experi- 
ments (radiative Feshbach spectroscopy) . The last exper- 
iments have been instrumental in observing I — 4 g-wave 
Feshbach resonances with ^mG resonance width. 

Based on the previous work [l5j . we have also im- 
proved the model for the Cs-Cs collision and in addi- 
tion used multi-channel bound state calculations to as- 
sign each Feshbach resonance with pertinent quantum 
numbers. The quantum numbers correlate each reso- 
nance to a molecular bound state at zero magnetic field. 
We identify 3 s-wave, 6 p— wave, 32 d— wave, 10 /—wave 
and 10 5— wave Feshbach resonances and one shape reso- 
nance. One resonance in the (3, 3) -I- (3, 3) channel could 
not be identified and is possibly a three-body collision 
resonance or a two-body Feshbach resonance with very 
high partial-waves number. 

The model has been used to calculate the molecular 
energy structure below the {F = 3) + {F = 3) contin- 
uum and the collision properties in the (3, 3) -I- (3, 3) and 
(3, —3) -I- (3, —3) scattering channels. These data will pro- 
vide crucial information for experiments on cesium Bose- 
Einstein condensation js^l and cesium molecules 0- 

In general, this paper presents results of a successful 
collaboration of experimental and theoretical efforts and 
resolves the collision anomalies of cesium atoms. The 
excellent agreement on >60 resonances between exper- 
iment and theory with only 5 parameters marks a tri- 
umph of the predictive power of atomic interaction the- 
ory. Since the experimental determination of the reso- 
nance locations is better than the theoretical estimate. 
Further improvements would require more flexibility in 
the (short-range shape) of the two Born-Oppenheimer 
potentials. Moreover, theoretical modeling will need to 
use thermally averaged elastic and inelastic rates as well 
as an improved model of the radiative lineshape. 
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